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Abstract 

We describe Lemaitre-Tolman-Bondi cosmological models where an anisotropic pressures is considered. 
By using recent astronomical observations coming from supernova of la types we constraint the values of 
the parameters that characterize our models. 
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I. INTRODUCTION 



About a decade ago, current local measurements of redshift and luminosity-distance relations of 
Type la Supernovae (SNIa)[l| indicate that the expansion of the universe presents an accelerated 
phase . In fact, the astronomical measurements showed that the SNe at a redshift of z ~ 0.5 
were systematically fainted, which could be attributed to an acceleration of the universe. These 
local observations can then be extrapolated to the universe at large, and the resulting Friedmann- 
Robertson- Walker (FRW) model can be used to interpret all cosmological measurements, which 
strongly support a dominant dark energy component, responsible for this cosmological acceleration. 

It is usual to say that this accelerated expansion is due to a cosmological constant, the so called 
concordance ACDM, which fits the observations quite well, but there is no theoretical understanding 
of the origin of this cosmological constant or its magnitude (the well known cosmological constant 
problem). Also, this model presents a coincidence problem of late cosmic acceleration, which 
establishes the fact that the density values of dark matter and dark energy are of the same order 
precisely today [J]. 

In order to get rid of the above problems there has been put forward other kind of different 
models, where some exotic, unknown and uncluttered matter component, dubbed dark energy^] 



(see also Refs 



IB, 



.a 



k-essence 



for recent reviews 



10h12||. phantom field 



is considered. Among these possibilities are quintessence 
13-15], holographic dark energy 16-18], etc. (see Ref. 
for model-independent description of the properties of the dark energy and Ref. [20|] for possible 
alternatives). However, although fundamental for our understanding of the evolution of the uni- 
verse, the nature of this quintessence remains a completely open question nowadays, moreover, 
very often they require fine-tuning 
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All these models work under the assumption of the Cosmological Principle, which imposes 
homogeneity and isotropy. But, the question is, are these symmetries consistent with observations? 
We know that inhomogeneities are abundant in the universe: there are not only clusters of galaxies 
but also large voids. It has usually been argued that these symmetries should only be valid on very 
large scales. 

It is true that this recognition is a result of the strict use of a FRW homogeneous and isotropic 
model. However, inhomogeneous models have been put forward. For instance, if phantom energy 
is considered, then it is possible to sustain traversable wormholes in which the phantom matter 



is an inhomogeneous and anisotropic fluid 



2314251] . Also, the authors of Ref. 26J] considered in- 



homogeneous metrics, but with an equation of state given by P = oop, with oj the equation of 



2 



state parameter being constant. Here, they study the gravitational collapse of spherical symmetric 
physical systems. On the same line of reasoning, the author of Ref. [27] studied a Lematre-Tolman- 
Bondi (LTB) model, taking into account a perfect fluid, in which they describe either a stellar or 
galactic system, where the collapse of these systems might yield to a possible formation of a black 
hole. Also, the effect of pressure gradients on cosmological observations by deriving the luminosity 
distance redshift relations in spherically symmetric, inhomogeneous space-times endowed with a 
perfect fluid was considered in Ref. 28]. 



If inhomogeneities are properly considered it might be possible to explain the observations 
without introducing the dark energy component. Recently, this possibility has renewed interest 
in inhomogeneous cosmological models, especially in the LTB solution 2914311] . which represents a 
spherically symmetric dust-filled universe. Because of its simplicity this solution has been consid- 
ered to be most useful to evaluate the effect of inhomogeneities in the observable universe, like the 



luminosity distance-redshift relation 
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■ 3j] . In most of these kind of models the matter component 
has been considered to be a perfect fluid with vanishing equation of state parameter, i.e. u = 0, 
corresponding to a dust fluid. 

Recently, numerical simulations of large scale structure evolution in an inhomogeneous LTB 



model of the Universe was studied in Ref. 



351 ] and the LTB model whose distance-redshift relation 



agrees with that of the concordance ACDM model in the whole redshift domain and which is 
well approximated by the Einstein-de Sitter universe at and before decoupling and the matching 
peak positions in the CMB spectrum was considered in Ref. 3a ] . Also, the kinematic Sunyaev- 



Zeldovich (kSZ) effect and the CMB anisotropy observed in the rest frame of clusters of galaxies 
was developed in Ref. [37]. 

In this work we want to address the study of an inhomogeneous LTB model, in which the matter 
component will be an inhomogeneous and anisotropic fluid 23j, |2j], where the radial and tangential 
pressures are given by p r = top and p± = u:±p, respectively. As far as we know no one has studied 
this problem previously. 

The outline of the paper is as follows. The next section presents a short review of the Lemaitre- 
Tolman-Bondi model. In the Sections III and IV we discuss the analytical solutions for our model, 
respectively. Here, we give explicit expressions for the functions of the radial coordinate and the 
luminosity distance. Finally, our conclusions are presented in Section V. We chose units so that 

c = h = i. 
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II. THE LEMAiTRE-TOLMAN-BONDI MODELS 



Lemaitre-Tolman-Bondi metric can be written as Refs. [29143 1| 



ds 2 = dt 2 _ ^Ytf dr 2 _ A{ t) 2 dU 2 (1) 

1 — k(r) 

where A is a function of the radial coordinate r and the time coordinate t, i.e., A(r,t), k(r) is an 
arbitrary function of the coordinate r, A'(r, t) = dA/dr and dVl? = d6 2 + sin 2 Odcf) 2 . 

Let us suppose that the universe is filled with fluid with a stress energy tensor as T^ u = 
diag(p,p±,p\\,pu), where p(r,t), p± and p» are respectively the energy density, the radial pressure 
P±( r it) = Pr(r,t) and lateral pressure pn(r,t) = p,p(r,t) = Pe(r,t) as measured by observers who 
always remain at rest at constant r, and 9. 

From the Einstein equations we obtain for the 00-component, we get 

A A' Mr) A 2 k'(r) 
A A' A 2 A 2 AN 

The other components of the Einstein equations; the rr equation gives 

A Mr) A 2 

2 A + A^ + A 2 ~ = - 87TGp± > (3) 
and from the 66 and (jxfi equations are 

^-(2AA + k(r) + A 2 )=-8nG Pll ^(A 2 ), (4) 

where the dot denoting the partial derivative with respect to t. 
Combining both equations ([3]) and ([H) we obtain 

From the conservation equation T fl u . fl = 0, we have that 

^ + ^(p± + p) + 2H( P + Pll ) = 0, (6) 

where H = A/ A. Note that if we restrict ourself to the case A(r, t) = a(t)r, these equations are 
exactly those obtained in Ref . . in which ^ + H{2>p + p±_ + 2p\\) = 0, with the identifications 
Pr = Pi and pi = Note also that when p± = pn, Eq.([5]) result in that p± is only function of 
time, and leads to the usual conservation equation. 

Now we shall require that the radial and the lateral pressures have barotropic equations of state. 
Thus, we can write p± = to±p and = wyp, where co± and wy are constant eq. of state parameters. 



Note that in the case when u± = cju = 0, we obtained the standard LTB model where the matter 
source have a negligible pressure, representing dust. 
Introducing the function 



F(r,t) = A(A 2 + k{r)) 

equations ((2j> and ([3]) can be written in a compact form 

F' = 8irGpA'A 2 , 
F = -8ttGuj ± pAA 2 . 



(7) 



(8) 
(9) 



The function expressed by Eq. ([7]) is the generalization of the one used in [43J , [44| and [45] , where 
the dust version were studied. We can recover this case by explicitly taking u± = in the second 
equation, which implies that F is only a function of r. 

In general, equation ([7]) could be considered as a Friedmann Equation for a non-zero pressure 
case. Thus, we rewrite 



H 2 (r,t)=H 2 (r) 



Vm(t)[^ + (1 - fl w (r)) ( 



(10) 



where we have defined as before H = A/A and 7 is a parameter that is zero for the dust case. The 
matter density and spatial curvature are defined by 

F (r) 



n M (r) - 1 



k(r) 
420-2' 



(11) 
(12) 



where the subscripts correspond to present values, Fq = F(r, to), Hq = H(r, to) and Aq = A(r, to). 
These relations are very important because, once we perform the analysis of our models with the 
observational data, we can actually find the profile fij^(r) that best fist the data. This is done 
explicitly in the next section. 

For an observer located at r = 0, incoming light travels along radial null geodesies, i.e. ds 2 = 
dQ 2 = 0, so time decreases with, dt/dr < 0, and thus we have 



dt 
dr 



A'(r,t) 
yJl-Kr)' 



(13) 



which together with the redshift equation (see Ref.[46|) 



dr 



y/l-k{r) 



dz (l + z)A'(r,tY 



(14) 



enable us to write down the relation between redshift z and time 

dt A'(r, t) 



(15) 



dz (l + z)A'(r,t) 

Equations (I14p and (I15D are used to find the functions t(z) and r{z). From these expressions, 
one can immediately write down the luminosity distance, the co-moving distance and the angular 
diameter distance function of redshift [45] 



d L (z) = (l + z) 2 A(r(z),t(z)), (16) 
dc(z) = (l + z)A(r(z),t(z)), (17) 
d A (z) = A(r(z),t(z)), (18) 

where z and A are evaluated along the radially-inward moving light ray. 

III. SEPARATION OF VARIABLE 

In the following we will derive analytic solutions for our LTB model, from separation of variable 
of the p and A, respectively. From Eq.(|5]) we find that 

p(r,t)=f(t)A- 2{1 -^\ (19) 

where f(t) is an arbitrary function of t. This suggest us that we can look for solutions under the 
following assumption: p = p\(t)p2{r), and also A(r,t) = A\{t)Ai{r). 
Inserting these identities in to Eqs.© and (0) we find that 

pi(t) = Ci^7 a , a = 3 + 2cj|| (20) 
P2(r) = C 2 A^, p = 2 (l-^y (21) 

Now, from Eqs.([2]) and ([3]) we find the following differential equation for the arbitrary function 
k(r) as a function of ^b( r ) 

2uJ± dAj + Aj <yl+ ^ = ~ Al ^ + 3UJ± ^ ~ 2AlAl = °' ^ 

As we can see, the left hand side depends exclusively of functions depending on the variable r and 
the right hand side everything depends on t, so both sides are equal to a constant C. We will have 
different family of solutions depending on the value that the constant C could take. Notice also 
that u)u is not relevant parameter in obtaining expressions for A\ and A 2 . 
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Using this separation of variables, the metric ([I]) takes the form 

~ dA 2 (r) 



ds 2 = dt 2 - A((t) 



A z 2 (r)dW 



(23) 



1 - k(r) 

where is clear that A 2 {r) can be considered as a new radial coordinate. In this case, instead of trying 
to get the explicit functional form r(z), we just need to obtain A 2 (z). Notice that although Eq. (|23l) 
looks very similar to the FRW one, the freedom in choosing k(r) keep the model inhomogeneous. 
However, if we assume in Eq. (|22p that uij_ = 0, then we get that k = CA 2 , and thus it reduces 
Eq. (|23p to the FRW homogeneous metric solution. It means that this separation of variables, 
reduces directly to the FRW model in the case of the universe it filled by dust. 

At this point, we can immediately write down the explicit expression for A±, in terms of the 
redshift z. From Eq. (fT5|) . we get 

dz (l + z)A 1 V ' 1 + z' V ' 

with Ai(fS) = A\{z = 0). Actually, this result is independent of any time dependance that we can 
anticipate for A\(t). This is somehow the analog to the FRW relationship between the scale factor 
a(t) and the redshift z. 

On the other hand, in general terms it is not possible to write down an explicit expression for 
A 2 (z). However, we could write down Eq. ([T4"]) under the assumptions set of the previous section, 
leading to the equation 

A' 2 dr = rflog(l + z) 

v 7 ! - W) A 

In order to solve this equation, we need the explicit expressions for Ai(t), and a relationship 
between k(r) and A 2 (r). Both relations follows from Eq. (|22p . In the next section, we consider 
several cases where this procedure enable us to write down an explicit expression for A 2 (z), which 
together with Eq. (|2"4"|) enable us to compute the luminosity distance given by Eq. (fT6j) . and thus, 
we could test each model by using supernova data. 

In the following we will study different families of solutions, which are characterized by the choice 
of the constant C (see Eq. ()22p ). In each case, we perform a Bayesian analysis using supernova 
data to constraint the parameters appearing in the model. If the corresponding scenario is favored 
by the observations, then we can use the results given by Eqs. fjlip and (|12p . and plot the profile of 
the function Qm (f). 

By using the separation of variables described at the beginning of this section, we found that 
H(r, t) = A\jA\ is a function of time only, so that Hq is a constant independent of the r coordinate. 



Also, if we take into account Eq. (|24p . we find that A$(r) = A\[z = 0)A 2 (r), where again it is clear 
that A 2 can be considered as a new radial coordinate. In the following, we proceed in order to find 
some specific solutions to our model. 



IV. SOME SPECIFIC SOLUTIONS 



The simplest case to solve Eq. (l22j) is when the constant C takes the value zero. Here, we obtain 
that 



Ai(t) = Cl [3(l + L0 ± )t}^+"±\ (26) 
A 2 (r) = 02^)'^!, (27) 

where c\ and c 2 are two arbitrary constants. The solution expressed by Eq. (|26p is obtained imposing 
the requirement that A\(t = 0) = 0. Combining Eqs. (|24p and (|26p we can find a relation between 
the time t and the redshift z given by 

where to expressed the present times. 

On the other hand, from Eq.([27]) we get that k(A 2 ) = c^2 (1+w±)/w± , and from Eq.([T2]), we 
obtain the following expression for the density parameter 

Note that A 2 could be considered as a new radial coordinate. This latter expression, allows us to 
plot the radial profile of the inhomogeneous density matter. Notice also that a "void" solution has 
to satisfy that as A2 00, Q,m has to increase towards an asymptotic constant value. This means 
that the second term in Eq. ([29|) have to decrease as A 2 — > 00. This is possible only if u± < — 1/3 
or uj± > 0. 

In order, to obtain A 2 (z) we need to integrate Eq. (|25p . In this way, the integral in the left hand 
side can be written as 

A * dr - f dA * (30) 

where a = — (1 + <^j_)/w_|_. This latter integral can be done analytically only in some cases where 
the parameter oj±, takes some specific values. 
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Case 



Here, we have a = 1 and the integral of the Eq. (l30p can be done directly. Inserting this 
expression in to Eq, (|25j) . and having in mind the Eqs. (124p and (|26l) . we get 



A 2 (z) = - 
c 



1 



l-cA 2 (0) 



4^1(1 + ^)4 - l}V 
2Ai(0) l V ' ij 



(31) 



The luminosity distance given by Eq. (|16p is then computed using the Eqs. (|24p and (|3ip . This 
expression has three free parameters; e\ = Ai(0)/c, e 2 = Cj4 2 (0) and = 3io/4. To test our model 
with observational data we use the Supernova Cosmology Constitution sample [^j], consisting of 



48J, 



397 SNIa expanded in a redshift range 0.015 < z < 1.5 and the more recent Union2 sample 
consisting of 557 SNIa. 

In Fig. [I] we see the confidence contours for the three parameters under consideration using the 
Constitution set. The best fit value for our parameter Ai(0)/c is shown with two horizontal lines 
indicating the range for 68.97% (continuous line) and 90% (dashed line) confidence region. 
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FIG. 1: Here, it is shown the result of the bayesian analysis of the model. We use the Constitution ( 47]) 
data set consisting in 397 SNIa. The best fit value for our parameters e\ — Ai(0)/c, e2 = c^2(0) and = to 
is shown with confidence contours showing the 68.97% and 90% regions. 



Using the numerical value for the parameters and the relation (|12|) we can find an explicit 
expression for the density parameter as a function of A 2 . Clearly the solution with lo± = — 1/2 
satisfy the requirement for a void solution, as we mention in the last paragraph. The results of 
the numerical analysis gives a very small value for ^(O), which does not permit us to estimate 
the value for Qm( z = 0). The same happens using the Union 2 data set [48|]. The results of the 
fitting process is shown in Fig.©. In the first case Xred = 1-18 and in the second one xf-ed = 0-98, 
showing that this is a quite good model. 
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FIG. 2: The figure shows the result of the bayesian analysis of the model. We use the Union 2 ( 48j) data 
set. The best fit value for our parameters e\ — A±(0)/c, e2 = 0^2(0) and = to is shown with confidence 
contours showing the 68.97% and 90% regions. 



B. Case oj 



Here a = 2 and from Eq, (|30l) . we get 

1 



M{z) 



sin 



2 Cl 



log(l + z) 



1 



sm 



Veto 

Mo) 



log(l + z) 



(32) 



Here, the solutions will be a sin or sinhfunction, depending on the sign of the c constant. Although 
it is possible to write down in explicit form the luminosity distance, and then test it under a bayesian 
analysis, we get that k (A2) oc A\, and thus leading to a FRW metric form and no inhomogeneity 
is obtained, as can be seeing from Eq, (|29p . Thus, in this case we recover the homogeneous FRW 
model. 



C. Case uj±_ = — |. 

Another analytical solution can be obtained by using ojj_ = — |. In this case a = 1/2 and the 
relation k{A2) oc VA2 enable us to have a inhomogeneous universe. This also satisfies the criteria 
of a "void" solutions, because uj± < —1/3. However, the bayesian analysis indicate that this model 
is not appropriated for describing the observations of the SNIa. 

Let us consider now the case where the equations expressed by ([22]) are equals to a constant C 
different from zero. The general solution for A2{r) it is found to be 

'+"1 p a2 



where the first term corresponds to the homogeneous solution of Eq. (|22p . i.e., C = (see Eq, (j27 
By taking into account Eqs. ([33|) and (fT2|) we notice that this class of solution are not permitted, 
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due to the presence of the second term proportional to A\ in Eg. (f33j) . We see that, as A2 — > 00 
the matter density parameter Qm goes to ±00, depending on the global sign of this term (which 
depends on C and u±). This behavior makes the model unacceptable. 



V. CONCLUSIONS 

In this paper we have studied a certain class of inhomogeneous LTB universes, in which the 

Tin 

matter component is an inhomogeneous and anisotropic fluid [23, 24], where the radial and tangen- 
tial pressures are given by p r = up and p± = uj±p, respectively. We have found analytical solutions 

to the system in order to compute explicitly the luminosity distance that was tested with the SNIa 

r I 

observations. We use both the Constitution 47] and Union 2 [48[] data sets. Our equations enable 
us to find directly, without any arbitrary ansatz, the profile of the density matter as a function 
of the radial coordinate. In this way, the fitting process with observations can fix the value of 
the density matter, both at z = and at z — > 00. We have found that the only model favored 
by the observations is of the "void" type solution, where VLm is a function that increase with the 
increment of the radial coordinate, leading to a constant asymptotic value. Although we can not 
determine the current value for the matter density, because of the very low value for ^.2(0) in the 
fitting process (see comments at the end of subsection IIV|) . we think that the method could be of 
interest in the case of nearly FRW inhomogeneous models. 

A necessary next step in testing these kind of models is to consider other observational probes. 
At present there is no agreement about how to proceed in the use of CMB and Large scale structure 



as additional probes 



38, 



39| even in the simplest LTB dust model. Although progress has been 
made [40|, we need to understand better the evolution of density perturbations in LTB spacetimes. 
Because at present we do not have a clear framework to address these problems, we have to look 
for additional observational probes. For example, in a recent work [41] the authors used galaxy 
ages [42]. The results still favored the ACDM model against the simplest dust LTB solution, in 
agreement with other probes. 

Note also that we have not addressed the stability issue of our models. But, since it becomes 
similar to a FRW- type of model (see expression (j23j) ) we expect that it model studied here is 
effectively stable. The demonstration could be carried in a similar way than that followed in the 
FRW case. However, here the situation become more complicate due to its homogeneity. We 
expect to address these points in the coming future. 
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